Modifications to the pulsar kick velocity due to magnetic interactions in dense plasma 
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In this work we have studied the non-Fermi liquid (NFL) behavior that enters into the expression 
of the pulsar kick velocity of the neutron star composed of degenerate quark matter core. We have 
incorporated leading order (LO) and next to leading order (NLO) corrections to the velocity and 
compared the results with the Fermi liquid case. Results for large magnetic field case persistent 
in neutron stars are also presented. The relation between radius and temperature has been shown 
with different values of the kick velocity with and without NFL corrections. 



Souvik Priyam 




* souvikpriyam.adhya@saha.ac. in 




2 



I. INTRODUCTION 



Exploration of the phenomenon of pulsar kicks i.e. the observed large escape velocities of neutron stars (NS) out 
of supernova remnants has drawn significant attention in recent years [l|-l5(- Understanding the origin of such high 
velocities, the so called kick velocities and its connection with pulsar spin has been the subject of many theoretical 
models. These include postnatal electromagnetic boost and hydrodynamical instabilities in the collapsed supernova 
to mention a few 0, 0). However, one of the most natural explanations for these large escape velocities is due to 
anisotropic neutrino emission from quark matter present at the core of the NS. It is well known that neutrinos carry 
away almost all of the energy released during the supernova and following proto-neutron star evolution. Thus an 
asymmetry of roughly about 3 percent in neutrino momentum will be sufficient enough to generate such phenomenal 
velocities |8[ . For a NS with quark matter core in normal phase, the main mechanism of neutrino emission has been 
the quark direct and inverse URCA processes given by [H, HJ , 



The above two reactions, in presence of a strong magnetic field can give rise to asymmetric neutrino emission as 
explained in [H, [l(J. The underlying mechanism for this is the possible polarisation of the electrons leading to 
neutrino emission in a preferred direction. For this to happen, the magnetic field strength has to be equal to or 
greater than some critical value which is given as B cr i C = mfc 3 /(qih), where m,i and qi are the mass and charge of 
the corresponding particle. Thus, we can easily compute that B & cric ~ 4.4 x 10 13 G and B q cric ~ 10 2 x B e cric G i.e. the 
value of the critical magnetic field for electron and quark system respectively. Such strong magnetic field forces the 
electrons to occupy the lowest Landau level polarising the electron spin opposite to the direction of the magnetic field. 
This leads to the formation of the neutrino and anti-neutrino emission cones which give rise to polarised neutrino 
emission along the axis opposite to the magnetic field direction. 

The electron polarisation for different conditions of magnetic field and kick velocities has been studied recently 
by Sagert et. al.Q. In this work, the authors have discussed pulsar acceleration mechanism based on asymmetric 
neutrino emission from quark direct URCA process in the core of the NS. Here the dependence of the kick velocity 
with the quark phase temperature and radius of the star with varying quark chemical potentials has been studied. 

One of the interesting devolopement in recent years has been the study of the non-Fermi liquid (NFL) behavior, 
of quantities like neutrino emissivity [O, [HI, mean free path (MFP)[ll|, Ell or the specific heat 0,0 of the dense 
quark matter which significantly modifies the initial cooling rate of the NS via neutrino emission. The NFL effects are 
dominant in the rclativistic regime driven by transverse or magnetic interactions which in the non-relativistic limit 
is p — v I 'c suppressed fbH ]. Infact the quark dispersion relation in dense plasma changes significantly in presence of 
magnetic interactions mediated by exchange of transverse gluons [HI, [l7| ■ It is also known , that the quasiparticle 
decay width which is given by the imaginary part of the quark self-energy also receives leading order (LO) corrections 
via the transverse interactions while the electric interactions mediated by the longitudinal gluon exchange contributes 
to the higher order. Furthermore, the damping rate (7) in such a plasma also goes as 7 ~ (E — fx), where E is the 
energy of the quasiparticle and \x is the corresponding chemical potential. It is to be noted that in the non-relativistic 
plasma, 7 ~ [E — /1) 2 [lril. fl7|. The quark self energy has been calculated recently both for the case of zero and low 
temperature by involving corrections up to next to leading order (NLO) terms which modifies the quark dispersion 
relation appreciably (l8| . Accordingly some of the physical quantities, like heat capacity, pressure are also modified 
[l5j |. In many calculations, such corrections appear in the modification of the phase space giving rise to appreciable 
contributions to the well known Fermi liquid (FL) terms. Such corrections are already been seen to be important 
resulting in enhanced neutrino emission from dense quark matter (Til fl2l | . The cooling behavior of NS with quark 
matter has also been studied recently by several authors H, [HI [H|. In all these cases, the correction at the LO has 
been seen to involve Tln(l/T) term which has been dubbed as anomalous corrections in many of the recent literatures 
[Til H7I . [HI . Similar behavior has also been reported while investigating quantities like drag , diffusion coefficients 
and thermal relaxation time of electrons for relativistic degenerate plasmas (l9l - [2l| . 

It is in this context, we revisit the problem of the calculation of the kick velocity to see whether such NFL corrections 
are significant to alter the kick velocity compared to the FL results. 

The plan of the paper is as follows. In section II we discuss the formalism to calculate the kick velocity and calculate 
the velocity for different conditions of the magnetic field. We present the results in section III followed by conclusion 
in section IV. 
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II. FORMALISM 



Neutron star is composed mainly of neutrons with a fraction of protons. Proton fraction may depend on various 
models. In presence of magnetic field, the proton fraction may increase [13]. We, however focus on the core of the 
NS where the density is expected to be so high that the relevant degrees of freedom may be considered to be quarks 
rather than hadrons. The pulsar kick velocity can be calculated if the luminosity, radius of the NS and the neutrino 
cmissivity is known. The emissivity is related, as we shall see, with the specific heat of the core. Initially, total 
neutrino luminosity (L) can be approximated as Q, 



L 



(3) 



where R is the radius of the NS and e is the neutrino cmissivity. Thus, due to momentum conservation we have, 
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In the above expression, the polarisation fraction of the electrons has been denoted by \ an d mass of the NS by Mjvs- 
Therefore the kick velocity can be written as H , 
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Using the cooling equation, 



C v dT = —edt 
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one can rewrite Eq(5]in terms of the specific heat of the quark matter core. In the region of low temperature and high 
chemical potential, the Fermi liquid expression of the specific heat of quark is given as [l4[ , 
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where N c and Nf are the number of color and flavor factors respectively Thus the Fermi liquid contribution to the 
pulsar kick velocity as reported in Q can be recast into the following form, 
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The calculation of such kick velocity performed so far is restricted to the Fermi liquid results [8[ . Currently, however 
it is shown in a series of works that NFL behavior arising out of magnetic interactions in the relativistic regime can 
contribute significantly and wins over the electric or Coulombic interaction for the case of ultradegenerate matter as 
mentioned in the introduction. In general, in a medium the quark dispersion relation can be obtained by solving the 
following equation [17j . 



w = (E p{uj) + ReS(w,p(o;))) 



(9) 



where the £ represents one loop quark quasiparticle self energy. Explicitly for excitation near the Fermi surface this 
has been derived in [171 UM , 
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Here e = to — fi q ~ T; u> being the quasiparticle energy; m 2 = (N f g 2 fj, 2 ,) / '(47T 2 ) and is related to the Debye screening 
mass by m 2 — m 2 D /2. In the above calculation of the quark self ener gy, the dominant contribution at one loop 



order arises when the gluon in the loop is soft (~ g/i) as shown in FigfT] 23 
gluon propagator to incorporate Debye screening and Landau damping 16l 



This approach requires dressing of the 
1?||. This modification to the in- medium 
dispersion relation results in a modified phase space which is responsible for appreciable corrections in all the cases 
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FIG. 1. Quark self-energy with resummed gluon propagator. 



including the specific heat of the degenerate quark core. Thus the specific heat of the degenerate quark matter up to 
NLO is given by [UGH, 
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where the coupling constant g is related to <? e // as, 
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and C, : 



total 



is the sum of the FL, LO and NLO contribution to the specific heat of the quark matter. Thus we obtain 



the LO contribution to the kick velocity as, 
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where Cf = (N% — l)/(2N c ) and the constants are c\ = —0.13807 and c 2 = 0.0530516. Now we have also extended 
our calculation beyond the LO in NFL correction. The NLO correction to the pulsar kick velocity is obtained as, 
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where the constants are evaluated as, 
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The net contribution to the pulsar kick velocity upto NLO is obtained by the sum of the Fermi liquid result and the 
non-Fermi liquid correction upto NLO: 
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The long range magnetic interactions lead to an anomalous r 2 lnT _1 term in the expression of the pulsar kick velocity. 
The specific heat of the quark matter is modified due to strong magnetic field present at the core of the NS. However 
it has been found that the effect of such high magnetic field is insignificant compared to the NFL corrections to the 
specific heat of the degenerate quark matter core in this temperature regime. 



A. Pulsar kick for the case of small magnetic field 



In presence of magnetic field (B) the number density of electrons is given by: 
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The number of Landau levels is limited to rj max = ^ ^ e . For the case of (n 2 — m 2 ) 3> 2eB, the number of occupied 
Landau levels is large. Thus we obtain, for the case of vanishing temperature (cold neutron stars), the electron spin 
polarisation is Q , 
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where the critical value of the magnetic field is given by B e cr ~ 4.4 x 10 13 G. Thus including the effect of electron spin 
polarisation on the FL contribution to the kick velocity, Eq.© is modified as, 
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The LO contribution to the kick velocity is given by, 
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We have also extended our calculation beyond the LO in NFL correction. The NLO correction to the pulsar kick 
velocity turns out to be, 
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The total velocity is given by Eq. fTT))) . 



B. Pulsar kick for the case of large magnetic field 

In this section the magnetic field strength is chosen to be much larger than the temperature, the chemical potential 
as well as the electron mass (/i e , m e , T <C \/2eB). The number density in this case is given by, 

The electron polarisation is given as @ , 
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Thus, the pulsar kick velocity is obtained as, 
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The LO contribution to the pulsar kick velocity is obtained as 

= 4N c N f ^ 2 q R 3 
V lo 27 \M NS 

2 5/4 



(C F a s ) 

T 3/2 ln 



Tin 



0.04 5 ^ 9V /iV7 
f 



dT 



ln(2) y VeB 

The NLO contribution to the kick velocity is obtained as 
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The above equations have to be solved numerically. The numerical estimation has been presented in fig. [5] . The total 
kick velocity up to NLO is given as the sum of the FL, LO and NLO contribution to the kick velocity. 



III. RESULTS AND DISCUSSIONS 



An estimation of the radius of the NS as a function of temperature (R-T relationship) has been presented in this 
section for different kick velocities and different conditions of magnetic field in NS. For this purpose, we have assumed 
the quark chemical potential to be AOOMeV with temperature ranging from 0.1 to 5 MeV. Our parameters are in 
good agreement with the high baryon density prevalent at the core of the NS. In addition, we have taken a s = 0.5 and 
/i e = 10 MeV Jj|. The pulsar mass has been taken to be 1.4M (where M Q = 2 x 10 30 fcg is the mass of the Sun). As 
high magnetic fields are persistent in the core of the NS, we have taken magnetic fields of ~ 10 16 G or higher relevant 
to the corresponding cases. In writing Eq.®, we purposely ignore the effect of magnetic field on C„ @. It has been 
seen that the effect of the magnetic field on C v is rather small compared to the result without the magnetic field. 

In the left panel of Fig[2] we observe that there is a considerable decrease in the radius of the NS with the inclusion 
of the LO correction over the FL case. This entails a significant increment of the kick velocity for a given temperature 
due to NFL correction. We have obtained the result considering fully polarised electrons, thus making the electron 
polarisation fraction (x) equal to unity. We have plotted the R-T behavior for two values of the kick velocities i.e. 
lOOOfcm/s and 500km/ s. In the right panel of FigJ5]we have presented the R-T comparison for the case of partially 
polarised electrons. In this case, the electron polarisation fraction \ is l ess than unity. 

In Fig J3] the left panel shows the comparison between the FL, LO and NLO result for the radius and temperature 
dependence for the case of fully polarized electrons for kick velocities of lOOOfcm/s. The right panel shows similar 
comparison for kick velocity of 500km/ s. 

The left panel of FigQ] shows the comparison between the FL, LO and NLO result for the radius and temperature 
dependence for the case of partially polarized electrons for kick velocities of 1000km/ s. The right panel shows similar 
comparison for kick velocity of 500km/ s. The LO correction lowers the radius of the NS for a specific value of 
temperature. We have found out that the NLO results impose a slight correction to the LO results in both FigO and 
Figg] 

Fig (5] shows for a large value of magnetic field, we have generated the plot numerically as this case is too complicated 
to be solved analytically. We have found that there is an appreciable change of the R — T relationship of the NS due 
to LO. However, the NLO correction to the kick velocity does not show considerable change from the LO result as 
evident from the graphs. 



IV. CONCLUSION 



In this work, we have derived the expressions of the pulsar kick velocity including the NFL corrections to the specific 
heat of the degenerate quark matter core. The contributions from the electron polarization (x) f° r different cases has 
also been taken into account to calculate the velocities. In addition, comparison has been made between the NFL 
LO and NLO contributions to the kick velocity with the FL case. We have found that the NFL LO contributions 
are significant while calculating the radius-temperature relationship as seen from the graphs presented for the case 
of the neutron star with moderate and high magnetic field. The anomalous corrections introduced to the pulsar kick 




FIG. 2. The left panel shows the comparison between the FL and NFL LO result for the radius and temperature dependence 
for the case of fully polarised electrons. The right panel shows similar comparison for partially polarised elctrons. The graphs 
have been plotted for kick velocities of lOOOfcm/s and 500km/ s for both the cases. 
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FIG. 3. The left panel shows the comparison between the FL, NFL LO and NFL NLO result for the radius and temperature 
dependence for the case of fully polarized electrons for kick velocities of lOOOfcm/s. The right panel shows similar comparison 
for kick velocity of 500km/ s. 



velocity due to the NFL (LO) behavior increases appreciably the kick velocity for a particular value of radius and 
temperature. However, for all the cases, no appreciable change in the R-T relationship has been observed for the 
NLO correction with respect to the LO case. 
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FIG. 4. The left panel shows the comparison between the FL, NFL LO and NFL NLO result for the radius and temperature 
dependence for the case of partially polarized electrons for kick velocities of lOOOfcm/s. The right panel shows similar comparison 
for kick velocity of 500km/ s. 
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FIG. 5. The left panel shows the numerical comparison where high magnetic field has been taken into account along with 
vanishing temperature for kick velocity of 1000fcra/s. The right panel shows similar comparison for kick velocity of 500km/s. 
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